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■ Abstract 

, , We develop a general low-energy multiple-scattering partial-wave theory for gapless topological 

C^ ' 

"^ . insulator (TI) surfaces in the presence of magnetic impurities. As applications, we discuss the 

c/3 ■ 

D I differential cross section (CS) dA/dip, the total CS Atot, the Hall component of resistivity fi, and 



inverse momentum relaxation time Tm for single- and two-centered magnetic scattering. We show 
that differing from the nonmagnetic impurity scattering, s— wave approximation is not advisable 



'^ ' and convergent in the present case. The symmetry of CS is reduced and the backscattering occurs 

O ' and becomes stronger with increasing the effective magnetic moment M of single magnetic impurity. 

We show a non-zero perpendicular resistivity component fi, which may be useful for tuning the 
^ , Hall voltage of the sample. Consistent with the analysis of dA/d(p, by comparing Tm with Ajot, 

o 

we can determine different weights of backscattering and forward scattering. Similar to CS, fi and 
'^ ■ Tm also exhibit oscillating behavior for multiple magnetic scattering centers due to interference 

rn 

O . effect. 

(N 



PACS numbers: 72.10.-d, 72.10.Fk, 73.20.-r, 73.50.Bk 



*zhang_ping@iapcni. ac.cn 
^ sslee@senii . ac. en 



I. INTRODUCTION 

A topic of fundamental importance in condensed matter physics is how the presence 
of defects or impurities induce strong modifications on the local electronic properties of 
crystalline solids. These modifications, with the stunning development of scanning tunneling 
microscopy (STM), have been extensively investigated on metal surfaces where they are 
well known as Friedel oscillations and manifest as standing waves in the local electronic 
density spanning regions up to ~10 nm from the defects on the metal surfaces. One kind of 
particularly suitable prototype that have been used in a large amount of STM measurements 
to study the effects of impurity and the formation of adsorbate superstructure are the (111) 
surfaces of noble metals, on which the surface-state electrons form a two-dimensional (2D) 
nearly free-electron gas. These Shockley-type surface states are dispersed as e = tt^k"^ /2meff 
(measured relative to the bottom of the surface-state band) and localized in narrow band 
gaps in the center of the first Brillouin zone of the (lll)-projected bulk band structure. 



Thus they have extremely small Fermi wave vectors kf = ^P2Jm^JJep j ft and consequently 
the Friedel oscillations of the surface state have a significantly larger wavelength than those 
of the bulk states. 

Recently, topological insulator (TI) has attracted tremendous experimental and theoret- 
ical studies [l|, |2|. Unlike (111) surfaces of noble metals, a peculiar characteristic of TI is 
the presence of strong spin-orbit coupling (SOC), which results in a variety of unique prop- 
erties. One intriguing fact is that the ideal TI surface is described at low energies by a 2D 
massless-Dirac wave equation with an additional locking between momentum and spin of 
surface electron. Because of the Dirac spectrum and SOC induced fermionic chirality, the 
impurity scattering effect in TIs is naturally expected to display novel behavior that should 
be absent from the conventional semiconductor or metal-surface 2D electron gases. Many 
experimental and theoretical efforts towards this issue have been payed. The anomalous 
Friedel oscillations in the vicinity of a single localized impurity (sMo], as well as the identifi- 



cation of the nature and the precise location of impurities on TI surface using STM [THuJ, 
have been discussed. However, when impurities are located close to each other, multiple 
scattering effects should be important, such as the issue of the long-range interactions be- 



tween the adsorbates mediated by the Dirac electrons of TI surface 12|, [iSj . In particular. 



since the quasiparticle's spin is strongly coupled to its momentum, quantum interference be- 



tween different spin states during multiple scattering process could then display prominent 



phenomena such as electric conductance weak (anti-) localization ij] and Aharonov-Bohm 



effect [15| in STM signals. 

In the presence of the time-reversal symmetry (TRS), the backscattering induced by 
nonmagnetic impurities is forbidden on the gapless TI surface because of a vr Berry phase 
associated with the 2tt adiabatic rotation of Dirac electron spin along the Fermi energy 
surface. However, considering the magnetic impurities on the gapless TI surface, one would 
like to observe the backscattering since the TRS is broken. Many efforts have been devoted 
to exploring this issue. For example, very recently, quasiparticle interference induced by a 



magnetic Co adatom on gapless Bi2Se3 surface has been found in STM experiments |16 |. 
Furthermore, a magnetic field can be generated when the TI sample is deposited on a 
lithographically patterned ferromagnetic layer, which could also induce backscattering of 
massless Dirac electrons in TI 17 1. 



Because of its importance both from basic point of interest and to Tl-based chemical 
catalysis and electronics, in the present paper we address this issue by presenting a first 
attempt at the theoretical evaluation of the multiple scattering problem of the massless 
Dirac electrons on the TI surface in the presence of localized and identical magnetic impu- 
rities. Specially, we present the analytical expressions for multiple partial-wave scattering 
of massless Dirac electrons with magnetic impurities, based on which the asymptotic mul- 
tiple scattering amplitude for random arraying magnetic impurities are obtained under the 
particular large distance approximations (the identical impurities are treated as a large 
scattering center). The differential and total cross sections (CSs), the inverse momentum 
relaxation time, and the Hall component of resistivity are discussed. We find that differ- 
ing from the nonmagnetic scattering, for the magnetic impurity scattering, the CS is not 
convergent under s— wave approximation. Therefore, higher partial waves should be intro- 
duced. For the single magnetic impurity scattering, we show the fact that the backscattering 
becomes much stronger when increasing the effective magnetic moment M. By comparing 
the inverse momentum relaxation time with total CS, we can determine whether there exist 
more backscattering than forward scattering or not. Similar to CS, the inverse momentum 
relaxation time and Hall factor display oscillating behavior for multiple magnetic scattering 
centers due to interference. 



II. MODEL AND THEORY 



The eigenstates of the effective low-energy Hamiltonian of TI surface near the Dirac point 

[18| 

Ho=hvf{crxk)-z (1) 



are given by the spinors \l/o,± (f) —~^ ( 1, =F^e*^'= ) , where the upper /lower sign corre- 
sponds to the electron/hole part of the spectrum. Here, f / ~ 5 x 10^ m/s is the Fermi 
velocity, cr are Pauli matrices, and k is the in-plane wavevector. The impurity potential can 
be expressed as 

V,=^J,Si-CTe{a-r), (2) 

where Si=Sni is the classical spin (with its orientation vector rii) of the ith magnetic 
impurity, Jj is the exchange coupling strengths, a is the radius of the scatterer, and 9 (r) 
is the Heaviside function. If the measurements are performed at a temperature higher than 
the Kondo temperature 191], the coupling between impurity spins will not exceed the critical 
value Jcr before a Kondo effect occurs. In this work we assume that the exchange coupling 
Ji < Jcr, so that the Ruderman-Kittel-Kasuya-Yosida interactions between impurity spins 
and Kondo screening of the impurity spin by the band electrons are neglected, and the 
impurity spin acts as a classical local magnetic moment under mean-field approximation 

m 

In order to obtain the analytical expressions of wavefunctions, we just consider the com- 
ponent of the classical spin along the normal line of TI surface, i.e., Vi=MazQ (a — r). The 
fact that a magnetic Co impurity with only perpendicular spin component on the TI surface 
does not open a gap has been experimentally observed [16|. To develop a scattering theory 
from localized, cylindrically-symmetric scatterers, it is convenient to resolve the problem in 
cylindrical coordinates. By considering the continuity of the wavefunction at the boundary 
of the magnetic scattering potential, one can immediately obtain the analytical expression 
of the scattered wave, written as 

oo 

^sc (r, k, ±) =soG+oTo<^± + J2 [siG+iT-+s-iG-iT^^] <^'^. (3) 

Here, $!|? denotes the incident plane- wave centered about a single scatterer located at r„. 



The cylindrically-symnietric Green's functions take the form 









2 \^ ^^W (^^^) g..„ ^(1) (^^^) g..„ J ' 

for |e| > M, where upper and lower signs in the right side of these expressions denote the 
e > and e < parts of the spectrum, k = ^, p^=r—rn and e^^" = ^" ^^ . For the energy 
regin of |e| < M, the Hankel functions H^ (kpn) in G±i should be replaced by the modified 
Bessel functions of first kind Ii (kpn). The scattering amplitude is expressed as 

A+Ji {k'a) Ji^i (ka) — A_Ji {ka) Ji^i {k'a) 
A.hI^^ (ka) Ji+i (k'a) - A+HJl\ (ka) Ji {k'c 

for |e| > M, where A± = a/Ic ± M\, k' = -'^—j- , and Ji is the Bessel function of order /. 

Whereas s; should also be replaced by 

Aj^Ii {k'a) Ji+i (ka) + A_Ji (ka) Ii+i {k'a) 
-A_hI^^ (ka) h+i {k'a) - A+hII\ (ka) h (k'a) 

for the case of for |e| < M. Note that si (si) satisfies the unitarity condition Re[s;] = — \si\ 
(Re[si] = - \sif), and 

lim si = lim si = - Jz+i (ka) /Hi+i (ka) (8) 

for all /. The /^''-partial-wave t-matrix is T,=^=diag( Pj^^ TiPl^i ) wi^^ Pj^ = ^^{dr±lde)'' ■ 
A detailed derivation is given in Appendix A. 

It is easy to extend scattering theory of massless Dirac fermions to the realistic and 
reasonable case of multiple magnetic impurities, where the quantum interference effect in 






the propagation process of Dirac fermions on T 



discussed in previous studies, such as Ref. [l7|. Taking into account all of partial waves. 



surface can be observed. This has not been 



for N magnetic scatterers located at positions T'i,r2, ■ ■ ■Tat, the scattered wavefunction is 
given by 

^,c(r-) = G(r)SD-V. (9) 

Here, G (r) (a 2x2A^ (2/jnax+l) matrix) contains the propagation information from detector 
to impurities G±i (r, rj). S* is a diagonal matrix with nonzero element s^/. The G matrix. 



is constructed by T^" [G±ii (n, m)], describing the propagation between impurities. can 
be written as a 2N (2/max+l) xl vector which imposes informations of incident waves (see 
details in Appendix B). 

At this stage, we should point out that the above equations for multiple magnetic scatter- 
ing of massless Dirac quasiparticles are similar to those for multiple nonmagnetic scattering 
of massive Dirac quasiparticles, but totally different from those for multiple nonmagnetic 
scattering of massless Dirac quasiparticles, since the expressions therein can be simplified as 
a more compact form due to s'_(i+i) = S;- 

The above theory enables to solve multiple magnetic scattering problems in gapless TI 
surfaces with higher partial waves, which could be important as distances between scatterers 
decrease or the scattering potential is strengthened. One simple application is to calculate 
the magnetic scattering CSs. To calculate the CSs, we have to take the approximations 

_ ^iiip -j?Qj, large distance. As 



Pn 



r-{x+iy) 



H (z) -^ J^e'y'-^-^), and e 
a result, 

vI/,e(r)^/(fc,V^)4^| , I, (10) 

V2r y^ie'^P J 

where / {k, ip) is the scattering amplitude, from which we have the differential and total CSs 
as follows: 

^ = l/(fc..)r. (11) 

A„ = r d^\f(k.^f= sj^lm [e-""f (k, ,3 = 0)]. (12) 

Here, we have used the two-dimensional optical theorem. 

Besides, we could also obtain the transverse component of resistivity (or say the analog 
of Hall component in the case with magnetic field) 

^= dip\f {k,ip)f sin ip, (13) 

Jo 

and the inverse electron momentum relaxation time (the quantity proportional to the dissi- 

pative component of resistivity) 

'.27r 



M 

'0 



[ d^\f{k,p)\'{l-cosp). (14) 

Jo 



III. RESULTS AND DISCUSSIONS 



In the following calculations, without losing the general properties, we shall just consider 



the incident wave $! 



V2 



1, -I 



propagating along the positive x direction. In par- 



ticular, for a single magnetic impurity scattering, we can obtain the scattering amplitude 
including all of the partial waves, which is written as 



/ (fc, ^) = /o (^) + Y. [/^ (^) + /-' (^)] e 



i(k—kf)-r' 



(15) 



1=1 



with /o (v?) = 


y ink^o and/±i(^) 


by 
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d(p irk 
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-s-t^e^*''^. The differential and total CSs are given 



So 



1=1 



sie 



ilip 



s-ie 



-il(p~\ 



Re (so) + /^ Re {si + s_ 



«=i 



(16) 



(17) 



This total CS equation is obtained from the optical theorem. It is clear that the s— wave is 
independent on the direction of scattered wave, therefore, we have to introduce higher partial 
waves, such as p— and rf— waves, and so on. If one just considers the s— wave in calculations, 
differential CS may lead to an unreasonable result of dA/dip\ ^^ = (backscattering is 
forbidden) for some particular effective magnetic moment M. This is different from the 
nonmagnetic impurity scattering on TI surface [12] as well as on conventional 2DEG with 
weak Rashba SOC [20|], where the s— wave approximation should be a reasonable choice. 

The results of normalized differential CS as a function of energy e for the massless Dirac 
electron scattered by a single magnetic impurity absorbed on TI surface with effective mag- 
netic moment M=60 meV are shown in Fig. [T](a). In the calculations we take /max = 2, 
which works out convergent results. Different from the nonmagnetic impurity scattering 
case, the backscattering is obvious in the differential CS, i.e., dA/dip\ ^^ ^ 0, in present 
case since the time-reversal symmetry is broken by magnetic impurity scattering. Further- 
more, we find that for the weak effective magnetic moment M (such as the values chosen in 
this work M < 100 meV), the backscattering is greater than forward scattering. However, 
if the effective magnetic moment is large enough (for example, when M ~ 500 meV and 
e = 450 meV) we find the backscattering is weaker than the forward scattering (not shown 
here) . 
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FIG. 1: (Color online) The normalized differential CSs dA/d(j) for (a) a single and (b) two magnetic 
scatterers located at ri^2 = (i3, 0) on TI surface with effective magnetic moment M = 60 meV. 
The radius of scatterer a=l nm, and /max = 2 are chosen. 

Two-impurity scattering provides a good test-bed to highlight the coexistence of various 
scattering phenomena, including transmission, reflection, interference, and resonance. The 
corresponding CS offers a measure of interaction events between the two impurity centers, 
and interference effects are useful in revealing actual electron density currents on TI surfaces. 
For instance, if two impurities are close to each other, the electronic wavefunctions will be 
scattered from both impurities, resulting in quantum interference. From the above theory, 
we can obtain the simple expression of scattering amplitude just containing the s— wave 
(^max = 0), which is given by 






n,m=l 



\/2mk 

\l J(2n-l),(2m-l) "*" L J 2n,2n 



(18) 



± 



[fl-] 



2n,(2m-l) 



[D-] 



(2n-l),2m 



However, if higher partial waves (/max > 1) are taken into account, the scattering amplitude 
expression for N magnetic impurities becomes tedious and complex since D— matix becomes 
a lager one. 

Typical numerical results of differential CSs for two magnetic impurities located at ri 2 = 
(±3,0) are presented in Fig. [l](b). We also note that the interference effect is related not 
only to the effective magnetic moment M but also to the configuration of impurities. For 
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FIG. 2: (Color online) The magnetic moment M dependence of differential CS dA/d(j) along the 
negative direction of x axis (/> = vr for one magnetic scatterer (red line), two scatterers (green 
line) located at ri^2 = (±3,0), and three scatterers (blue) located at ri = (—3,0), r2 = (3,-1), 
r^ = (2,4). The Fermi energy is chosen as e = 70 meV. 

the present considered impurity locations ri 2, on one hand, we find from Fig. [Hb) that the 
backscattering is more prominent than the forward scattering. On the other hand, comparing 
with the nonmagnetic double-impurity scattering on gapless TI surface, the symmetry of 
differential CSs for two identical magnetic scatterers is reduced. 

On one hand, independent on the impurity locations, with increasing the effective mag- 
netic moment M, we find that the relative strength of backscattering becomes more and more 
remarkable since the differential CS along the negative x direction dA/d(p\ ^^ increases with 
M, see Fig. [2] with the Fermi energy e = 70 meV. On the other hand, comparing with the 
scattering from a single (red curve) magnetic impurity, two (green curve) or three (blue 
curve) impurities will weaken or strengthen the backscattering due to interference effect, 
which is dependent on the impurity configurations relative to the direction of incident wave. 
We must point out that in the calculations, we should use s/ for the energy region of |e| < M, 
while Si for the energy region of |e| > M, which are denoted in Fig. |2l 

Now let us turn to discuss the total CSs, which are exhibited in Fig. [3l As mentioned 
above, the s— wave approximation cannot give out convergent result, see the black curves 
in Fig. [3l whereas, when we introduce higher partial waves (such as /max = 2 chosen in our 
calculations), the total CSs becomes convergent ultimately. Differing from the nonmagnetic 
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FIG. 3: (Color online) The total CS Ktot for a single (a-b) and two (c-d) magnetic impurities 
located at ri^2 = (i3, 0). The effective magnetic moment is chosen as M = 60 meV in (a)and (c), 
while M = 100 meV in (b) and (d), respectively. 

impurity scattering, although the higher partial waves can induce remarkable corrections, we 
have not found additional resonant peaks in total CSs due to higher partial waves. Moreover, 
it is obvious that interference between double impurities brings about oscillations in total 
CSs (see Figs. Ws) andEl^d)), which cannot be observed in the case of nonmagnetic impurity 
scattering 12[. From numerical calculations, we find on one hand that, the optical theorem 
is correct and should characterize the general multiple-scattering processes, since the results 
obtained by the numerical integration of the first equality in Eq. flT2|) are in good agreement 
with that obtained from the second equality; On the other hand, the curves of total CSs are 
smooth at the energy of e = M, which indicates that the limiting function of si and s/, i.e., 
Eq. (IH]) is reasonable. Besides, for the much strong effective magnetic impurity scattering, 
with increasing the Fermi energy we find that the total CS for single- (double-) impurity is 
convergent to 4 nm (8 nm). This is also different from the nonmagnetic impurity scattering, 
where the total CSs converge to zero with increasing the energy of Dirac electrons [121]. 

In spite of the differential and total CSs, we also calculated the transverse component of 
resistivity fi, which is analogous to Hall component in the case with external magnetic field. 
The typical results of fi as a function of e for single- and double- impurity with different 
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FIG. 4: (Color online) Hall component of resistivity fi for a single (a) and two (b) magnetic 
impurities on TI surface; Inverse momentum relation time Tm (thick solid curves) and total CS 
Atot (thin dashed curved) for a single (c) and two (d) magnetic impurities as functions of energy e. 

M are listed in Figs. ID^a) and HJ^b), respectively. We find that the Hall component of 
the resistivity Q always keeps its sign as negative (i.e., fi < 0), which is independent on 
the impurity locations. In the numerical calculations, we take /max > 2 which results in 
a convergent result, however, taking into account higher partial waves, it is difficult to be 
obtained analytically from Eqs. (TT3|) and (TT5|) since the expression for integral result is 
tedious and complex. Therefore, the low-energy Dirac electrons are deflected to one side of 
TI sample due to magnetic impurity scattering. This fact may be helpful for tuning the Hall 
voltage of sample. Interestingly, we note that this type of Hall component also occurs in the 
nonmagnetic impurity scattering on gapless TI surface. We now switch gears and consider 
the case of double magnetic impurities on TI surface again. Different from the single case, Q 
exhibits oscillating behavior due to the interference during the multiple impurities scattering 
processes, as shown in Fig. HJ^b). 

Before ending this paper, we would like to discuss another important quantity, the inverse 
electron momentum relaxation time Tm, which is proportional to the dissipative component 
of resistivity. The numerical results are plotted in Figs. HJ^c) andlD^d) for single- and double- 
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impurity cases, respectively. The behavior of Tm (thick sohd curves) is quahtatively similar 
to the one for the total CSs Atot (thin dashed curves), and the interference effect is also 
clear in the Tm induced by double-impurity, see Fig. ID^d). The inverse electron momentum 
relaxation time Tm is a fairly sensitive quantity that determines whether the charge carrier 
is attracted to an impurity or is repelled from it. It is also useful for determining whether 
the backscattering is greater than the forward scattering by comparing it with the total CSs 
(note that both Atot and Tm have the dimension of an length in two dimensional scattering). 
If Atot<^M (Atot>rM), there is more (less) backscattering than forward scattering. Taking 
M = 60 meV as an example, we find Atot< Tm, as revealed by the red curves in Fig. lU^c), 
which indicates that the backscattering is greater than the forward scattering in the low- 
energy region. This fact is consistent with the behavior of differential CSs shown in Fig. 
[U^a). Particularly, by observing the green curves in Fig. lU^d) of energy region of 80~120 
meV (the shadow region), one can find Atot>^M- This suggests that due to the interference 
the backscattering is weaker than forward scattering in this energy region for the present 
double magnetic impurity locations, which is also consistent with the analysis of differential 
CSs. Consequently, we believe that our results shown here are reasonable, and we hope our 
findings could be detected in the future experiments. 

IV. CONCLUSIONS 

In summary, we have proposed a general low-energy multiple-scattering partial-wave the- 
ory for quasiparticles on the gapless topological insulator (TI) surfaces in the presence of 
magnetic impurities. Based on this theory, one can solve the scattering problems of A^ 
magnetic impurities. As an application, we have calculated the CSs, the inverse momen- 
tum relaxation time, and the transverse resistivity component for a single and two circular 
magnetic scattering. We have found that the usual s— wave approximation is not sufficient, 
while higher partial waves must be introduced to obtain convergent results. On the gapless 
TI surfaces, differing from the single nonmagnetic impurity case, the backscattering occurs 
and becomes much stronger with increasing the effective magnetic moment M. Interference 
effects are obvious in CSs from quasiparticle scattering off two magnetic scattering centers, 
and oscillating behaviors are introduced in Atot associated with higher-order partial-waves. 
A non-zero perpendicular resistivity component has also been shown. Similar to the total 



12 



CS, the inverse momentum relaxation time and the transverse resistivity component exihibit 
oscillations for multiple magnetic scattering centers due to interference. Furthermore, our 
theory could be extended to spin-polarized case. It could also be applied to simulate the 
electron flow (charge current and spin current) through a quantum point contact on a TI 
surface by monitoring the changes in conductance through the quantum point contact as a 
moveable STM tip is scanned above the surface of TI. 
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APPENDIX A: DERIVATION DETAILS OF SCATTERED WAVEFUNCTION 

Starting from the model considered in the main text, one can find the spinor spherical 
wavefunctions inside the magnetic scattering potential (r < a) as follows: 



.iw ( , /|e + M|iy/'''^ {k'r) 
^2 |e| k' \ ±^|e-M|i//|f {k'r) e'' 



Cr' (r, k', ±) = ^^== I -^^^^'^7, ,, :: '/ 1 , (H > M), (Al) 



we^'^ / J\^ + Mill (k'r) 

V2 |e| k' \ ^,/\^^M\Il+^ (k'r) e^' 



where k' = -^f- and w = e**^ is an overall phase meaning of the nonrelativistic wavefunc- 

tion in the rest frame of e = M. Here, HJ^ ' (z) and Ii (z) are /*^-order Hankel functions and 
modified Bessel functions of first kind, respectively. We would like to point out that since 
the wavefunctions should be zero at r = 0, we have neglected the modified Bessel functions 
of second kind Ki (z) for |e| < M, which are emanative in the limit of z — ?■ 0. Whereas, the 
wavefunctions outside the magnetic potential (r > a) can be expressed as 
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with k = -j^. x*^^^ (x*-^-*) denotes the outgoing (incoming) cyhndrical wave about r=0. The 
incident plane-wave centered about a single scatterer located at r„ is given by 



$- (r) = 5^ 



-^' (xS (Pn) +xS (Pj) 



(A4) 



where pn=r—rn and e^^" = ^" '^ ^^ . Remember that ^fc is the angle defining the direction of 
the wave vector, and 9^ is related to the direction of r—Vn in this equation. Then, the fully 
scattered wave function in the region pn > a can be written as explicitly 



^i(r,fc,±) = <l>i"(r) + ^sc(r,fc,± 

oo 

= y/ke'''-''" Yl 



l=—oo 



\e'''^xf\Pn.k.^) + \xf\Pn.k.^) 



ji(e„-ek) 



(A5) 



where 5i are phase shifts of the outgoing cylindrical partial waves, Xi {Pni^^^)- ^"^ ^^e 
region of p^ < a, the fully scattered wave function is given by 



j^/^ik-Vn y^°o 



i''di 



V(i) 



Cr^ (P„, fc', ±) - ICl'^ (p„, fc', ±)J e^'^^'-^'"), (|e| > M) 



^TT fr, fc', ±) 



(A6) 

By the continuity of the wavefunction at p„ = a, \l'i (a, fc, ±) = \l'ii (a, fc', ±), we can obtain 
the scattered wavefunction Eq. ([3]), and the scattering amplitude si and 'si shown in Eqs. 
(E]) and (^^ in main text. It is clear that the scattered wavefunctions shown in Eq. (jS]) are 
different from those for nonmagnetic impurity scattering on gapless TI surface, where 



*.c (r) = Y, 



°° l!^G,(r.r„,e)r;[4i"] 



1=0 



with 



WjWn 



Gi(x 



Hre 



Ti^Se-^C+i)^" 



±/7/!)e-'('+i)^" H^'^e-'' 



Wn 



(A7) 



(A8) 



w+i 



T;=diag( p-, P+ ), and 



Ji {k'o) Ji+i (ka) - Ji (ka) Ji+i {k'o) 



r(i) 



(1) 



Hl'^ (ka) Ji+i (K'a) - Hil[ (ka) Ji {k' 



(A9) 



Here k'=j-^, and Vq is the scalar potential. 
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APPENDIX B: DETAILS OF THE EXPEDITION FOR MULTIPLE SCAT- 
TERED WAVE 

In order to understand the extending operation, we would like to start from the s— wave 
scattering for two identical magnetic impurities. The total wavefunction can be written as 

2 

v]/ (r) = $ (r) + Y, 4"^G+o (r, r„, e) T,- [v[/„ (r„, k, ±)] , (Bl) 

ra=l 

where 

2 

^n (r) = $ (r) + 5^ 4™)G+o {r, r^, e) T,~ [v]/^ (r^, fc, ±)] . (B2) 

Equation (Bl) indicates that if the value of \E' (r) and its derivatives due to the /g ^ de- 
pendence of T^ at each scatterer is known, the entire wavefunction \1' (r) is completely 
determined. We can calculate the derivatives of \E'i at Vi and "$2 at r2 and combine the 
result into a matrix equation, which are given by 

T„-I*.(n)l\ ^_7r„-|*(n)l 



where 



with 



D — l4x4 — 64x484x4, (B4) 



To-[G'+o(r2,n,e)] 

s^'^ 

S= I ° (2) I ® 12X2. (B6) 

Sq 



Finally, the total wavefunction is written as 



^ (r) =$!"+G (r) SD-i ( ^° '"^ ^'''^^ I (B7) 

To [^ (r^)] 



with 



G(r-)=(G+o(r,ri), G+o (r, rs) ) • (B8) 

Taking into account /^ax > 1 partial waves, for A^ magnetic scatterers located at positions 
ri, r2, ■ ■ ■ Tat, the scattered wavefunction is given by Eq. ([9]) in main text. For numerical 
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calculations, we have to align the matrix elements reasonably, thereby, we define symbols 
/io = 2 (n - 1) (2/i„ax + 1) + 1, Ao = ;Uo + 1, ao = 2 (n - 1) (2/„,ax + 1) + 4/ - 1, 70 = ao + 1, 
To = 2 (n - 1) (2/,nax + 1) + 4/ + 1, r^o = To + 1, z/o = 2 (m - 1) (2/„,ax + 1) + 1, /3o = «^o + 1, 
a = 2 (m - 1) (2/,„ax + 1) + 4/' - 1, 7 = a + 1, z/ = 2 (m - 1) (2/,„ax + 1) + 4/' + 1, and 
P = u + 1. Following this way G (r) (a 2x2N (2/jnax+l) matrix) is aligned as 

(B9) 



G(r)= (G{r,r,), G{r,r2), ■■■, G'(r,rjv) 



with G (r, T-i) =[ G+o, G+i, G^i, ■ • ■ , G+z^,, G-^^, ]• Explicitly, for / = 0, we align 



'l.Mo 



^2,M0 



'l,An 



^2,Ao 



G+o{r,r„), 



(BIO) 



and for Z > 1, 



'!,«() 



'2,ao 



'l,n-, 



^2,ro 



'1.70 



Gfr: 



2,70 



'i,* 



^2,»7o 



The S* matrix is diagonal. 



S = diag ( s 



G+i{r,rn), 
G-i{r,rn). 



2A'(2«niax + l)x2Ar(2imax + l) 



with S, 






,(«) 



.(") 



(Bll) 
(B12) 



(B13) 



the G matrix, which is written as 



"no^m 



s_i . Then we construct 



/ 



G 



^(1,2) ■■■ G{1,N) 
G(2,l) ■■■ G{2,N) 

G{N,1) G{N,2) ■■■ 



[B14) 
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where G{n,m) is a 2 (2/inax + 1) x 2(2/max + l) matrix, which are constructed by 
Tj^ [G±i> {n,m)]. Exphcitly, one would ahgn the G matrix by the following way 



for / = and n ^ m, and 



G(/io,z/o) G(/io,/3o) 
G(Ao,z/o) G(Ao,/3o) 
G(/io,a) G(;Uo,7) 
G(Ao,a) G(Ao,7) 

G(/io, //) G(/io,/3) 
G(Ao,z/) G(Ao,/3) 

G{ao,uo) G{ao,(3o) 
G(7o,//o) G(7o,/3o) 

G(ao,a) G(ao,7) 
G(7o,a) G(7o,7) 

G(ao, '^) G(ao,/3) 
G(7o,//) G(7o,/3) 

G(ro,z/o) G(ro,/3o) 
G(?7o,t'o) G{r]o,l3o) 

G(ro,a) G(ro,7) 
G(?7o,a) G(r7o,7) 

G(ro,z/) G(ro,/3) 
G(r7o,^) G(r7o,/3) 



To" [G+o (n, m)] , 
^o~ [G+i' {n, m)] , 
^o" [G-i' {n, m)] , 



Tf [G+i' (n, m)] , 
Tf [G-i' (n, m)] , 
T;+ [G+o (n, m)] , 
T;+ [G+i' (n, m)] , 
T+ [Gt, (n, 



772 



(B15) 
(B16) 

(BIT) 



(B18) 
(B19) 
(B20) 
(B21) 
(B22) 
(B23) 



-> 



for / 7^ and n 7^ m, while G (n, n) = for n = m. (p can be written as a 2N (2/inax+l) x 1 
vector, 

^ ' H, 02, ■■■, 0^) , (B24) 



where 



To- [$ (r,)] , Tf [$ (r,)] , T+ [$ (r,)] , • • •, T,_ [$ (r,)] , T^_ [$ (r,)] 



Explicitly, for /' = 0, 



-> 



</> uo = Pq 



Ak 



•'•VV2 



-> 



^pr 



idkJk-Vn 1^/2 



Tie "e 



(B25) 
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For /' > 0, 



0a 


= Pv 


Q^k-r^l^ 


, 0^, 


= ^Pv+i 


^^g*e,g*fe.r„/^ 


(B26) 


(/^u 


= P^ 


gifc.r„/^ 


, 0/3 


= ^^/+1 


±ie'<^ke^k■r^/^ 


(B27) 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8] 

[9] 

[10] 

[11] 

[12] 
[13] 
[14] 

[15] 
[16] 



[17] 



M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011). 

Q. Liu, C.-X. Liu, C. Xu, X.-L. Qi, and S.-C. Zhang, Phys. Rev. Lett. 102, 156603 (2009). 

R. R. Biswas and A. V. Balatsky, Phys. Rev. B 81, 233405(R) (2010). 

R. R. Biswas and A. V. Balatsky, Phys. Rev. B 83, 075439 (2011). 

A. M. Black-Schaffer, and A. V. Balatsky, Phys. Rev. B 85, 121103(R) (2012). 

T. Zhang, P. Cheng, X. Chen, J.-F. Jia, X. Ma, K. He, L. Wang, H. Zhang, X. Dai, Z. Fang, 

X. Xie, and Q.-K. Xue, Phys. Rev. Leet. 103, 266803 (2009). 

P. Roushan, J. Seo, C. V. Parker, Y. S. Hor, D. Hsieh, D. Qian, A. Richardella, M. Z. Hasan, 

R. J. Cava, A. Yazdani, Nature 460, 1106 (2009). 

J. Wang, W. Li, P. Cheng, C. Song, T. Zhang, P. Deng, X. Chen, X. Ma, K. He, J.-F. Jia, 

Q.-K. Xue, and B.-F. Zhu, Phys. Rev. B 84, 2135447 (2011). 

Z. Alpichshev, J. G. Analytis, J.-H. Chu, I. R. Fisher, Y. L. Chen, Z. X. Shen, A. Fang, and 

A. Kapitulnik, Phys. Rev. Lett. 104, 016401 (2011). 

Z. Alpichshev, R. R. Biswas, A. V. Balatsky, J. G. Analytis, J.-H. Chu, I. R. Fisher, A. 

Kapitulnik, Phys. Rev. Lett. 108, 206402 (2012). 

Z.-G. Fu, P. Zhang, Z. Wang, and S.-S. Li, Appl. Phys. Lett. 99, 232109 (2011). 

Z.-G. Fu, P. Zhang, Z. Wang, and S.-S. Li, Phys. Rev. B, 84, 235438 (2011). 

H.-T. He, G. Wang, T. Zhang, I.-K. Sou, G. K. L Wong, J.-N. Wang, H.-Z. Lu, S.-Q. Shen, 

and F.-C. Zhang, Phys. Rev. Lett. 106, 166805 (2011). 

Z.-G. Fu, P. Zhang, and S.-S. Li, Appl. Phys. Lett. 99, 243110 (2011). 

M. Ye, S. V. Eremeev, K. Kuroda, E. E. Krasovskii, E. V. Chulkov, Y. Takeda, Y. Saitoh, K. 

Okanioto, S. Y. Zhu, K. Miyamoto, M. Arita, M. Nakatake, T. Okuda, Y. Ueda, K. Shimada, 

H. Namatame, M. Taniguchi, A. Kimura, phys. Rev. B 85, 205317 (2012). 

A. Zazunov, A. Kundu, A. Hiitten, and R. Egger, Phys. Rev. B 82, 155431 (2010). 



18 



[18] H.-J. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S-.C. Zhang, Nature Phys. 5, 438 

(2009). 
[19] J. J. Cha, J. R. Wihiams, D. Kong, S. Meister, H. Peng, A. J. Bestwick, P. GaUagher, D. 

Goldhaber-Gordon, and Y. Cui, Nano Lett. 10, 1076, (2010). 
[20] J. D. Walls, J. Huang, R. M. Westervelt, and E. J. Heller, Phys. Rev. B 73, 035325 (2006). 



19 



